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Abstract 

In this paper, connection between surface roughness and directed polymers in random medium are stud- 
ied, when the surface is considered as a directed line undergoing stochastic oscillations. This is performed 
by studying the influence of a stochastic elastic forcing term —Ky+rj(s), on aparticle moving along arough 
surface. Two models are proposed and analysed in this way: the random-walk process (RW) in its discrete 
and continuous form, and a Markovian process via the Ornstein-Uhlenbeck (O-U) process. It is shown that 
the continuous RW leads to an oscillator equation, via an effective action obeying a KPZ equation which 
is solved analytically. The O-U process allows to obtain information on the profile of surface for a long 
size substrate. The analogy with the roughness is achieved by introducing a quantity suited to directed line 
formalism: the height velocity variation dh/ds. 

PACS numbers: 05.40.-a, 02.50.Ey, 68.35.Ja 
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I. INTRODUCTION 



Characterization of surface roughness or interface growth is of first importance in many ap- 
plications: microelectronics, preparation of catalysts, magnetic materials (see e.g. yj]- IJ). And 
numerous models based on numerical as well as analytical approaches have been proposed over 
past decades to treat interface properties. These works concern in particular the KPZ equation 
(scale laws, roughness, growth) iS-l 15], and models treating connections between surface growth 



and directed polymers in random medium Ill6h - ll24ll . The present paper focus on the second class 



of these models, considering the roughness of a surface as the result of specific stochastic pro- 
cesses, and based on the directed line formalism. This formalism, much used in polymer physics 
[25|], allows to analyse the motion of a system by the way of transversal and curvilinear coordi- 
nates (26D-D28D (respectively noted y and s in this paper). And we propose models showing the 
analogies between the profile of a rough surface and the study of the transversal position equa- 
tion. Note that, in this formalism, s plays the role of the time. For example dy/ds designs the 
transversal velocity. Moreover, this formalism is also linked to the statistical physics: the partition 
function of the system and its links with the KPZ or the Burgers equation, which follow from this 
formalism, represent useful tool s Il29l]-[| 3lll. Concerning the roughness analysis, different works 
have treated this problem ||29ll . Il32ll - ll37ll . but the roughness is nevertheless poor studied from a 



purely stochastic point of view. Among recent works treating stochastic surface properties, we 
can cite [38] and 113911 which both have a Markovian approach, using a Langevin equation with a 
stochastic noise. 

In the models presented here, we assumed that the profile of surface is described by the tra- 
jectory of a particle undergoing a directing stochastic strength in addition to diffusion. We will 
use for this purpose the elastic forcing term f(y, s) = —ny + r/(s) inducing stochastic oscilla- 
tions, or equivalently, a potential such as / = —d&/dy; r] having for example the properties 
of a Langevin force. For deposition processes, we suppose that the surface start at the height h , 
and we are interested in the height velocity variation (hvv), dh/ds, which gives local information 
on the structure of the surface to analyse: for a fast variation of the height, we will have a high 
value of the hhv, (Fig. (Q])). Hence, this quantity is the reflect of local irregularities. But, the 
hvv do not have to be confused with the quantity dh/dx, where x stands for the horizontal coor- 
dinate: locally, the slope dh/dx may takes an infinite value (i.e parallel to the h axis), whereas 
{h(si) - h(s 2 ))/{si - s 2 )) is weak. 
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FIG. 1: Profile of surface: result of the conflict between the random medium and a directing strength 
—d^/dy. The local slope (a) (resp. (b)) stands for a weak (resp. high) value of the hvv. 

The outline of the paper will be as follow: we first derive properties related to profile of surface 
from pure and forced Random Walk (RW). In what follows, we will employed the term forced 
RW to design a RW undergoing a forcing term in addition to the diffusion. Then, we justify the 
use of an oscillator model comparing the directed line to a string vibrating. This will allow to use 
the restricted partition function to determine the transversal equation of motion, in presence of a 
quadratic stochastic potential. In the last part, we shall concentrate on the O-U equation of motion 
for the transversal position, and particularly on the asymptotic solutions for s — > oo. 



II. TOPOLOGY OF SURFACE VIEW AS A RANDOM WALK PROCESS 



The models presented here concern the use of the directed line formalism associated to a RW 
process. In its discrete form, we based our approach on the paper of Vilgis 12511 . treating directed 
polymers in random media. Our derivation is somewhat different of Vilgis, but leads to the well- 
known partition function of polymer physics as showed later on. 



A. Discrete random walk process 

We start from a set of vectors {bj = 6j e s ; i — 1, . . . , N}, which are successively connected, 
and recovering the whole of the surface (Fig. ©). We suppose they are statistically independent. 
The probability to find such a set is given by: 

P({b i ;i = l,... J /V}) = I]p(b i ), (D 

8=1 

where the probability p(bj) corresponds to a step of the random walk of bj, with the mean value 
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FIG. 2: Random Walk: N steps of length b are required in order to cover the surface on a linear size 
substrate L. 

properties: 

(hi) = (bib,-) = b%. (2) 

Sij is the Kronecker symbol and b is a characteristic length of the shift of position of the surface 
from one to the other. For the step bj, a normalized gaussian probability law is taken: 

p(b,) = (l/27r6 2 ) 1 / 2 exp(-b, 2 /26 2 ), (3) 

and the associated distribution of the surface length reads 

P(R) = Jdb 1 ...Jdb N 5(R-J2b^j Jlp{bi); (4) 

the vector R = J2iLi characterizing the line. Then, using the Fourier transform representation 
of the Dirac delta, a Gaussian law for the length distribution arises quite naturally 

P(R) = (l/27r6 2 iV) 1 / 2 exp(-R 2 /26 2 iV). (5) 

And from this expression, the quadratic mean "length" of the surface can be evaluated: 

(R 2 ) = j R 2 P(R 2 ) d R = b 2 N. (6) 

Note that the quantity (R 2 ) preserves information on the length of the trajectory, contrary to the 
geometrical quantity R. However, this result is a limiting case of a pure RW trajectory since it 
implies that the entanglements, and so the length of the surface, could escape to infinity. Thereby, 
for this case we label N such as 

(R 2 ) rw = b 2 N rw , (7) 
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with the meaning that iV rw steps are necessary in average for a pure RW to describe the profile of 
surface. The other limiting case, i.e. the straight line (si), is given by 

(R 2 ) sl = b 2 N 2 ~ L 2 , (8) 

in such a way that the general case (g), representing a forced RW, and included between the cases 
© and © , may be written 

(R 2 ) g = b 2 N^; N vw <N g < N sh (9) 

where 7 is a critical exponent representative of a strength acting on the RW. In addition, we have 
the conditions 

b 2 N 2 "<^b 2 N 2 for iV g — > N s \\ (10) 

which implies 7 = 1, and 

b 2 N 21 — > b 2 N rw for N g — iV rw ; (11) 

implying 7 = 1/2. Thus, 1/2 < 7 < 1, available for a line moving in the plane (h, x). This result 
may also be expressed in term of the length / ~ a/ (R 2 ) g of the surface: 

/ = b iV g 7 . (12) 

In the case of a two-dimensional RW, 7 may takes the form 7 = l/d{, where di designs the fractal 
dimension of the directed line. As a matter of fact, d f describes the degree of irregularity of the 
surface, and for 7 = 1 and 7 = 1/2, df becomes identical to the topological dimension of the 
system. Clearly, the exponent 7 is linked to the number of steps N g . A relation between the 
two quantities may be obtained from the probability of the general case, P g (R), when the line 
undergoes both effects of diffusion and of an external force. Anticipating the continuum approach 
described afterwards, we write this probability 

P g (R) = Afexp (-R 2 /26 2 iV g - /3$(R)) , (13) 

where $(R) stands for an external potential. For example, in the case of a harmonic potential, 
$ = (k/2)R 2 , the quadratic mean length is 
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which leads, with Eq. ©, to the expression 

ln(iV g /(l + 6 2 /3KiV g )) 
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(15) 



2 In (^ g ) 

One can then obtain information on the topology of a surface, in the sense that the ratio //Lis 
the reflection of its irregularities and roughness. The fractal dimension may be evaluated experi- 
mentally, for example, by the box-counting method with boxes of length b (in the same way as for 
electric discharges Q40D). In simulations, the varying parameter would be the number of steps of 
the RW. 



B. Equations of evolution 



If we express now the step as bj = Rj — Rj_i, and writing shortly P = P({hi, i = 1, . . . N}), 
the probability (Eq. (Q])) may be rewritten 



P 



1 \ N/2 



(2^) 



cxp 



1 N 

— ^2(Rt - Hi-iY 



2b 2 



8=1 



(16) 



One can then associated to the exponential argument a symbolic Hamiltonian 



[3H 



1 N 

— J^(Ri -Ri_i) 2 . 



2b 2 



(17) 



i=i 



This is the energy of a mass-spring chain of a one dimensional lattice. This suggest to take the 



continuous limit of this discrete chain: R ; — Ri 1 — > <9R/<9s, and YliLi ~~ * In > where s is the 



-N 



dimensionless curvilinear coordinate. In this limit, the probability of the pure RW reads 



AAexp 



2b 2 



N 



fdR\' 



yds J 



(18) 



where J\f is a normalization constant. From this result, one can construct the partition function 
as a sum over all the possible paths starting from s = and reaching the linear size substrate at 

s = N: 



Z = TV ex P 

All paths R(s) 



2b 2 



V 9s I 



(19) 



The argument of the exponential is equivalent to the transversal kinetic energy, or more precisely 
to the effective action of the system since we integrate over the "time" s. Furthermore, this path 
integral obeys the classical diffusion equation which arises naturally from a pure RW process: 



d_ 

ds 



Z(R,s) 



b 2 d 2 
2 dR 2 



Z(R,s), 



(20) 



with the boundary condition Z(R, 0) = 5(H). For a forced RW, we introduce a potential term 
$(R(s)) standing for an external (virtual) strength field, in such a way that the total effective 
action takes the form 



(3S 



2b 2 



N 



dR(s) 
ds 



ds + 



N 



$(R(s))ds. 



(21) 



We note S to emphasize the fact that this quantity represents a pseudo action (the integration is 
made over a Hamiltonian and not over a Lagrangian). This yields to the total partition function of 
the system, as the Feynman-Kac path integral: 



Z(R, s)=N [ exp - pS VR(s) 



(22) 



It is then convenient to introduce the transversal coordinate y = R/6 and so to write the action in 
the form |26] 



ds'H(y,s'). 



(23) 



2\ds'J 

H standing for the Hamiltonian of the system. By this way, we obtain the well-know expression 
of the partition function of directed polymers 



ry( s ) 

Z{y,s) =N /Pyexp 
Jy(o) 



-(5 / ds' 



ydyy 



v2VdsV 



(24) 



expression directly issued from a RW process with a directed line formalism. For a one dimen- 
sional system, Z(y, s) obeys the following PDE |4jJ]: 



§- s Z(y,s) 



v- 



d 2 Z(y, s) 



dy 2 



-P$(y,s)Z(y,s) 



(25) 



where v = 1/2(3 is a viscosity coefficient. 

Let us underline that there exists several versions of the Hamiltonian depending on system to 

study li^-lii, boh - n2.4n . 

In view of these results, several analogies may be drawn between directed lines and vibrating 
strings. In its discrete form, the directed line is composed of a succession of rigid chains and 
the whole evolving with a resulting motion depending on the form of the noisy potential. For 
the continuous form, the term of inertia in polymer physics is due to the thermal agitation of the 
medium, and the potential to an external field {e.g. electric field). In our case, the external potential 
is necessary because of the strength which have to be imposed to the line, in order to obtain realistic 
profiles of surface. So, the noisy potential of our model will be composed of a harmonic part and 
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of a stochastic part which shall make the oscillations stochastic: $(?/, s) = \ky 2 + r)(y, s), where 
k is the elastic constant of the line. The next step is to determine the equation of motion with the 
frequency of vibration of the line. In order to do this, let us proceed as follow: 

In one hand, starting from Eq. (l25l) . it is well-known that we may obtain successively a KPZ 
and a Burgers equation (approach (A)): 

OS d 2 S lfosY ^ 

and 

du du d 2 u <9$ ,„„ v 

— + u— = v— + — . (27) 

os dy ay* ay 

In this case the transformations are 

u = an d Z = exp (^S/2uJ . (28) 

Usually one introduce the free energy. But the free energy tends to the total energy of the system 
S when the entropy tends to zero. 

In the other hand, starting from a general Burgers equation for a velocity field u, we have 
(approach (B)) 

du du d 2 u dV 
as ay ay z ay 

and 

?L - v — + I (—\ 2 + v, (30) 

ds dy 2 2 \ dy J ' 



and finally, with Z = exp (F/2u), 



|-Z(x, S ) = z/V 2 Z(x, S ) + ^. (31) 
as 2v 

Whether u represents the transversal velocity, we have in this case 

u= T- = "IT- W 
as ay 

From these two approaches, a question arises: "May the quantity u of Eq. (|271) be interpreted as 
the transversal velocity of Eq. (|29l ) ?" In other words, is there identity between S and F ? 

To remove this ambiguity, let us start from the Burgers equation (|29l ), u given by (|32l ), and let 
us use the relation 

dF dF 

dF(y,s) = —ds + —dy; (33) 
8 



which yields to 



d 2 F 1 fdF\ 2 dF 
dF = v^ds + - — ds + V ds + — dy; 
dy 1 2\dy) ay 



and consequently to the expression 



(34) 



(35) 



The first and last integrals may be rewritten 



d 2 F ^ _ d_ f dF_ _ d_ f dy _ 
dy 2 dy J dy dy J ds 



and 



dF , 

dy 



dy, 
ds 



In such a way that the "free energy" will read 



F(y,s) = -v - 



1 /dy 



2 V ds' 
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ds' 



(36) 



(37) 



(38) 



Finally, the Hopf-Cole transformation Z = exp (F/2u), with l/2v = (3, yields explicitly to 



7 -V(y,«0 



(39) 



The heat equation (1251) being invariant under the shift Z ^ 2 = c x c 6 K, we wiil have 
equivalence between the two approaches if $ = —V. And, therefore, one obtain the classical 
Hamilton equation of motion: 



dy 
ds 



dF_ _dS_ 
dy dy 



fy 
ds 2 



dH 

dy 



(40) 



The determination of an equation of motion for y is henceforth equivalent to solve the KPZ equa- 
tion (1261) . Moreover, we have a relation between the height of the growing surface and the transver- 
sal position: \h(s) — h \ = |y(s) — y(0)|, and so, by extension, 



dy dh 
ds ds 



(41) 



For example dh/ds = dy/ds = 0. Thereby, the action may as well be written Q42 , 



43] 



S 



(42) 



We shall now concentrate on the solution of the KPZ equation (1261) . 



Analytic solutions of noisy KPZ equations have already been studied in the literature in specific 
cases I144ll - ll46ll . But the non-linearity of this equation added to the presence of a stochastic noise 
make the resolution difficult. Nevertheless, it is known that this equation may be transformed 
either into a Burgers equation or a Schrodinger equation with an imaginary time. And, among 
works which treat exact solutions of theses kind of equations, we can cite Il47ll - ll5 ill . This is our 
concern here to treat the case of a line undergoing the stochastic elastic forcing term — ny + r)(s), 
where the function rj contains all the noise that affects the oscillations (e.g. a Langevin force). This 
amounts to saying that we solve the KPZ equation (1261) with a potential $ given by s) = | y 2 — 
yrj(s). Moreover, we assume a constant initial condition for the transversal velocity: u(y, 0) = u . 

Following methods presented in Ref. 115 ill , it is shown in appendix that the solution is of the 
form: 

S(y,s) = ^-y 2 -y^(s) + Q(s), (43) 



where ip(s) and 0(s) are functions of s and of r/(s) (see the appendix). But, we have the 

relation ^ 

(hi <9? 

dJ = = "**(*) +#0- (44) 

Consequently, we obtain an equation standing for the equation of a forced oscillator with a term 
of friction: 

d 2 ll 

(45) 



d^ + m dS + uj{s)y = F ^ 



where the pulsation to 2 and the forcing term F are respectively given by d£/ds and dip/ds. We 
can notice from the expressions of £ and ip that, while the frequency / oc ^/d^/ds depends only 
on s, the term of force depends also on the stochastic noise 77. Thus, we expect profiles of the form 
drawn in Fig. ©. 

Once y determined, we can deduce the hvv from Eq. (|44l) : 



dh 
ds 



-£(s)exp f-J\(t)dt\ 



1 + / i){t) exp QT £(t')dt^J dt J + (46) 



Then, the degree of irregularity of a surface, and so its global roughness, may be estimated by 
adding up all the local slopes \dh/ds\, in such a way that the interesting quantity to determine is 
(\dh/ds\), where the mean is taken over all the length of the surface. 

Note that this problem is different of the one of a stochastic oscillator, which is a specific 



problem 11521 . 



53] 
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y 



(a) Harmonic 



(b) Stochastic: 
Weak Noise 



(c) Sochastic: 
Strong Noise 



FIG. 3: Variations of the transversal position y for a forcing term inducing oscillations depending on the 
intensity of the noise. Each mode (a), (b), and (c) leads to very different roughness. 

III. LONG SIZE SUBSTRATE 

We treat now the case of a particle of mass unity, moving along a rough surface when we study 
its profile at a fixed time. 

Let us consider the propagation of a Brownian particle through virtual scatterers. Its trajectory 
describing the surface to analyse. We assume that the particle obeys a general Langevin equation 
where an external force —d&/dy = —ky is considered: 

d^ = ^dS- Ky+m - (47) 

This means that the particle undergoes effects of the medium (term of friction) in addition with an 
oscillator term. Then, the Langevin force T(s) may be expressed in term of a Gaussian with noise 
r](s), such as 

= -C^ - k y + VWr,(s); (48) 

where V = £ 2 i>. By nature, the solution of this equation follows a Markovian process. When the 
inertial term is negligible, one obtain an Ornstein-Uhlenbeck process for the transversal position 

y- 

^ = -~y + V2fiTl(.8). (49) 



This well-known process 1541 15511 is associated to the following Fokker-planck equation for the 



transition probability P(y, s\y', s'): 



dP _ kdyP | _ d 2 P 
ds £ dy dy 2 
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With the standard initial condition P(y, s\z/q, sq) = S(y — y'), the solution writes [i56ll 

TTTTT^ ex P 



This probability is asymptotically "stationary": 

lim P(y,y',s) = 



k(y 



- exp 



2^(1 



&y 2 

2^v 



,-(3R/0«) 



(51) 



(52) 



In practice, the condition s ^> £/2k is sufficient to obtain this solution. This implies that, under 
this condition, the fluctuations of the transversal position y becomes independent of the curvilinear 
abscissa. Moreover, the quantity A = ^^v/k, appearing in Eq. (l52l) with the dimension of a 
length, takes the sense of the mean range for the transversal fluctuations with respect to a given 
position: 

r+oo 

(53) 



(54) 



y 2 P(y)dy = A 2 . 

> 

And the Einstein relation v = gives us the expression 

A oc k- l/2 



Thereby, the higher the elastic constant, the lower the range of fluctuations. The connection be- 
tween the O-U process and the roughness is then highlighted in Fig. (HI), where we can notice 
that, for s ^> 1, the maximal fluctuation values of the height h(y, s) are of the order of A. As a 
consequence, the hvv will be such as 

dh 



ds 



<A Vs>£/2£. 



Consequently, in this limit, 



dh 

ds 



)<A; 



(55) 



(56) 



available for models concerning processes of deposition executed from the left to the right of the 
figure, and on a "large" length of substrate. 

Then, an interesting point lies in a connexion between the O-U process and KPZ equations. 
Indeed, the transformation 

P(y,s) = Z(y,s)e-^ 4i> , (57) 
applied to Eq. (l50l) . leads to the heat equation 



dZ_ 

ds 



v 



d 2 Z 
dy 2 



k n 2 y 2 
2 ~ ~4zT 



(58) 
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FIG. 4: Asymptotic profile of surface for an O-U process. A has to be interpreted as the maximal range of 
fluctuations. 



and so to the KPZ equation 

dS „d 2 S 1 / dsY fk k 2 y 2 




ds dy 2 2 \ dy I V 2 4z> 



(59) 



And once more, the introduction of a quantity u such as u = —dS/dy, yields to a Burgers equa- 
tion: 

du ^ du _, d 2 u <9$ ^ 
9s <9y dy 2 dy 
But the same interrogation as the one of the RW process may be done, namely, may u be considered 

as the velocity dy/ds of Eq. (|49l ? As a matter of fact, these relations have been obtained with 

the assumption d 2 y/ds 2 ~ 0. So, let us assume that u = dy/ds. The consequence is that the first 

term of Eq. (|60l) vanishes, and so we obtain 

ldu 2 „d 2 u d§ 
v 



2 dy dy 2 dy 

^d_ fl /djA 2 + | > 1 9H _ fdy\ (61) 
dy I 2 \ ds / J dy dy 2 \ ds J 

Furthermore, we have that —dH/dy = d 2 y/ds 2 ~ 0. So the last part of Eq. (|6T| ), leads to a 
solution of the form 

^ = c 1 (s)y + c 2 (s). (62) 

Thus, u = dy/ds u obeys to Eq. (|62|) . Then, we refound Eq. (|49l by putting ci = — «/£ and 

c 2 (s) = v / 2^^(s). 

Once again, a stochastic process (here Markovian) may be connected to the coupled 
KPZ/Burgers equations, and by extension to the Hamiltonian of the system. In addition, a par- 
tition function like the one of the RW may be obtained with appropriated v and k. Moreover, this 
shown the relevant choice of a harmonic potential ~ y 2 to treat these kind of fluctuations, since a 
Markovian process and a RW process yields both to the same class of equations. 
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IV. CONCLUSION 



Two short stochastic models of directed line have been presented and applied to the description 
of rough surfaces. The background of our approach was that roughness of a one-dimensional 
surface is created by the motion of a virtual particle dived in a stochastic harmonic potential. It 
brings out of the continuous RW that we can rigorously linked the partition function of the line to 
the Hamilton equation of motion for the transversal position. The survey of roughness amounts 
therefore to solve a KPZ equation with a particular potential. Then, interpreting the roughness as a 
result of a stochastic elastic forcing term acting on a surface, we have obtained a forced oscillator 
equation where the force behaves as a noise for the line. The analytical expressions of the pulsation 
and of the force have been determined too. Consequently, local analysis of the roughness may be 
performed thanks to the expression of the hvv (|46|) which, besides, obeys a Burgers equation. 
For a sufficient long substrate, the O-U model shows that the maximal values of roughness are 
asymptotically bounded in a range of length A oc k~ 1 / 2 . The elastic constant being an adjustable 
parameter mastering the magnitude of the fluctuations. Moreover, we have shown that the two 
models, used with a same kind of force, had common properties since they lead both to the same 
class of equations. As a matter of fact, the two equations of motion are connected to the partition 
function of the system. Finally, these two models belong to stochastic processes described by the 
Burgers dynamics. 

Although we have presented one-dimensional models, these approaches could open interesting 
ways of investigation for numerical and analytical studies, in particular if one needs to generate 
surfaces with specific profiles or properties. 
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APPENDIX: SOLUTION OF THE KPZ EQUATION 



We solve here the KPZ equation (1261) or, equivalently, the heat equation (1251) written in the 
following way: 

d n d 2 Z ( (3k 2 



with the initial condition Z(y, 0) = exp (—u y/2u). The resolution is based on the Time Space 
Transformation method (TST) of 115 ill . We give only the main points of resolution. We would like 
to underline for the comprehension that the notations x and y are reversed between the present 



paper and H51H : y <-> x. 

The first step is to determine the initial condition. For this, we need to determine the expression 
for the variables a±, a2(s), p(s), q(s) and r(s) of [51]. We obtain successively: 

ai = \fk~/4v, (A.2) 

a 2 ( s ) = e-^ s (l + (3 rj^e^dtj , (A3) 

p( s ) = e ^ s } (A.4) 

q(s) = 2vs+ I [ rj^e^'dt'dt, (A.5) 
Jo Jo 

r = e^ s sinh( v ^s)/ v / ^. (A.6) 
In such a way that the resolution of (IA.1I) amounts to solve the equation: 

<9$ <9 2 $ 



dr dx 2 

with the initial condition 



(A.7) 



0) = exp (-aix 2 - (1 + u /2u)x) . (A.8) 

The solution reads explicitly 

s 1 io [ 9 , vr{2a^x + 1 + Ur ) l2v) 2 1 

$(x, r) = (4 ai UT + 1)~ 1/2 exp j-aix 2 - (1 + u /2u)x + \avT + 1) — J ' ( } 

Then, reversing the different transformations, we obtain as a solution of (IA.1I) 

s) = exp jy 2 f ai - aie 2 ^ s + y=e 2 ^ s tanh(v^s)^ + 
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y ^a 2 (s) - 2a x q{s)e^ s - (1 + u Q /2u)e^ s + ^^(1 + V 2 ^ + 2a 1 q( y s))e^ s tanh(v^s)^ + 
(^a iq 2 (s) + -^=(1 + 2a 1 g(s)) 2 tanh( v ^s) - - - In (cosh(^s))^ | (A.10) 



expression that we write for further analyse 



where 



f (s) = 4z/ f fll - a ie 2 ^ s + ^=e 2 ^ tanh(V^)) , (A.12) 
^( s ) = -2zy f a 2 ( s ) - 2aig(s)e^ s - (1 + M /2z/)e^ s + ^(1 + m /2^+ 



+2aig(s))e v/Ss tanh( v ^s)J , (A.13) 
e(s) = 2i/ ( -aig 2 (s) + -^=(1 + 2a!g(s)) 2 tanh(V«s) - - ^ In (coshes)) 

y y/ At A A 



(A.14) 

The quantity d, a 2 (s) and g(s) being respectively given by Eqs. (IA.2I) . (IA.3I) and (IA.5I) . Note that 
we recover the initial condition 0). Then, since Z = exp (^S/2v^ , we can write the solution 
on the final form 

S(y, s) = %V - w/>(s) + e(s) (A.15) 
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